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Abstract. A d- dimensional quantum model system confined to a general hypercubical geometry with linear 
spatial size L and "temporal size" 1/T ( T - temperature of the system) is considered in the spherical 
approximation under periodic boundary conditions. For a film geometry in different space dimensions i<7 < 
d < Icr , where < a < 2 is a parameter controlling the decay of the long-range interaction, the free energy 
and the Casimir amplitudes are given. We have proven that, if d = a, the Casimir amplitude of the model, 
characterizing the leading temperature corrections to its ground state, is A = — 16£(3)/[5cr(47r)' T / 2 .r(<7/2)]. 
The last implies that the universal constant c = 4/5 of the model remains the same for both short, as 
well as long-range interactions, if one takes the normalization factor for the Gaussian model to be such 
that 5=1. This is a generalization to the case of long-range interaction of the well known result due to 
Sachdev. That constant differs from the corresponding one characterizing the leading finite-size corrections 
at zero temperature which for d = cr = lisc = 0.606. 
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1 Introduction 

The confinement of quantum mechanical vacuum fluctua- 
tions of the electromagnetic field causes long-ranged forces 
between two conducting uncharged plates which is known 
as (quantum mechanical) Casimir effect [0,0]. The con- 
finement of critical fluctuations of an order parameter field 
induces long-ranged forces between the surfaces of the film 
HQ. This is known as "statistical- mechanical Casimir ef- 
fect" . The Casimir force in statistical-mechanical systems 
is characterized by the excess free energy due to the finite- 
size contributions to the free energy of the system. In the 
case of a film geometry L x oo 2 , and under given boundary 
conditions r imposed across the direction L, the Casimir 
force is defined as 



F, 



Casimir 



(T,L) = - 



8L ' 



where /° X (T, L) is the excess free energy 

f™(T,L) = f T {T,L)-Lf hx ^{T). 



(1.1) 



(1.2) 



Here f T (T, L) is the full free energy per unit area and per 
fcfiT, and /buik(r) is the corresponding bulk free energy 
density. 

The full free energy of a d-dimensional critical system 
in the form of a film with thickness L, area A, and bound- 
ary conditions a and b on the two surfaces, at the bulk 



critical point T c , has the asymptotic form 

f a , b {T c , L\d) Si Lf hulk (T c \d) + / s a m ,acc (T c \d) (1-3) 

+.fL iacc (T c \d) + L- {d - 1) A a , b (d) + --- 

as A — > oo, L> 1. Here /surface is the surface free energy 
contribution and A atb {d) is the amplitude of the Casimir 
interaction. The L d ependence of the Casimir term (the 
last one in Eq. ( |l.3| )) follows from the scale invariance 
of the free energy and has been derived by Fisher and de 
Gennes ||. The amplitude A a .^(d) is universal, depending 
on the bulk universality class and the universality classes 
of the boundary conditions . 

Equation (|1.3| ) is valid for both fluid and magnetic 
systems at criticality. Prominent examples are, e.g., one- 
component fluid at the liquid- vapour critical point, the 
binary fluid at the consolute point, and liquid 4 He at 
the A transition point jjj. The boundaries influence the 
system to a depth given by the bulk correlation length 
£oo(T) ~ |T— T C | _I/ , where v is its critical exponent. When 
Coo(r) <C L the Casimir force, as a fluctuation induced 
force between the plates, is negligible. The force becomes 
long-ranged when ^ (T) diverges near and below the bulk 
critical point T c in an O(n), n > 2 model system in the 
absence of an external magnetic field HQ]. Therefore in 
statistical-mechanical systems one can turn on and off the 
Casimir effect merely by changing, e.g., the temperature 
of the system. 

The temperature dependence of the Casimir force for 
two-dimensional systems has been investigated exactly only 
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on the example of Ising strips ||. In 0(n) models for 
T > T c the temperature dependence of the force has been 
considered in Q . The only example where it is investigated 
exactly as a function of both the temperature and mag- 
netic field scaling variables is that of the three-dimensional 
spherical model under periodic boundary conditions |^|. 
There exact results for the Casimir force between two walls 
with a finite separation in a L x oo 2 mean-spherical model 
have been derived. The force is consistent with an attrac- 
tion of the plates confining the system. 

The most of the results available at the moment are 
for the Casimir amplitudes. They are obtained for d = 2 
by using conformal-invariance methods for a large class of 
models Q]. For al ^ 2 results for the amplitudes are avail- 
able via field-theoretical renormalization group theory in 
4— e dimensions [§,[5||§], Migdal-Kadanoff real-space renor- 
malization group methods |10| , and, relatively recently, by 
Monte Carlo methods |tJJ| . In addition to the flat geome- 
tries recently some results about the Casimir amplitudes 
between spherical particles in a critical fluid have been 
derived too For the purposes of experimental veri- 

fication that type of geometry seems more perspective. 

It should be noted that in contrast with the quantum 
mechanical Casimir effect, that has been tested experi- 
mentally with high accuracy [ fl3| , the statistical- mechanical 
Casimir effect lacks so far a satisfactory experimental ver- 
ification (for comments on the specific difficulties that the 
experiment stacks with see, e.g. 

In recent years there has been a renewed interest 0, 
|l5[ in the theory of zero-temperature quantum phase tran- 
sitions. In contrast to temperature driven critical phenom- 
ena, these phase transitions occur at zero temperature as 
a function of some non-thermal control parameter, say g, 
(or a competition between different parameters describ- 
ing the basic interaction of the system), and the relevant 
fluctuations are of quantum rather than thermal nature. 
In the present article we consider a statistical-mechanical 
Casimir effect when critical quantum fluctuations play an 
essential role. 

It is well known from the theory of critical phenomena 
that for temperature driven phase transitions quantum ef- 
fects are unimportant near critical points with T c > 0. It 
could be expected, however, that at rather low (as com- 
pared to characteristic excitations in the system) tem- 
peratures, the leading T dependence of all observables is 
specified by the properties of the zero-temperature criti- 
cal point, say at g c . The dimensional crossover rule as- 
serts that the critical singularities with respect to g of a 
d-dimensional quantum system at T = and around g c 
are formally equivalent to those of a classical system with 
dimensionality d + z (z is the dynamical critical exponent) 
and critical temperature T c > 0. This makes it possible to 
investigate low-temperature effects (considering an effec- 
tive system with d infinite spatial and z finite temporal 
dimensions) in the framework of the theory of finite-size 
scaling (FSS). This theory has been applied to explore the 
low-temperature regime in quantum systems [ pT[^5| ,[l6|, 
when the properties of the thermodynamic observables in 



the finite-temperature quantum critical region have been 
the main focus of interest. 

In this paper a theory of the scaling properties of the 
free energy and Casimir amplitudes of a quantum spher- 
ical model [jlTj with nearest-neighbor and some special 
cases of long-range interactions (decreasing at long dis- 
tances r as \/r d+(T ) is presented. These interactions enter 
the exact expressions for the free energy only through their 
Fourier transform wich leading asymptotic is U(q) ~ q a , 
where a* = min(a, 2) [pi. As it was shown for bulk sys- 
tems by renormalization group arguments a > 2 corre- 
sponds to the case of finite (short) range interactions, i.e. 
the universality class then does not depend on a . Val- 
ues satisfying < a < 2 correspond to long-range in- 
teractions and the critical behaviour depends on a. On 
the above reasoning one usually considers the case a > 2 
as uninteresting for critical effects, even for the finite-size 
treatments . However recent Monte Carlo results sug- 
gest that it might well not be the case at least for contin- 
uous Ising model There Bayong and Diep state that 
for d = 2 the critical exponents does not depend on a and 
reach their short-range values for a > 3. On the basis of 
that result it seems that for finite-size systems the case 
a > 2 is nontrivial. Since, up to the authors knowledge 
that is the only example where a > 2 is of interest for 
studying critical properties, here we will consider only the 
case < g < 2. 

The investigation of the Casimir effect in a classical 
system with long-range interaction possesses some pecu- 
liarities in comparison with the short-range system. Due 
to the long-range character of the interaction there exists a 
natural attraction between the surfaces bounding the sys- 
tem. One easily can estimate that in the ordered state the 
L-dependent part of the excess free energy that is raised 
by the direct inter-particle (spin) interaction is of order of 
L~ a+1 . In the critical region one still has some effect stem- 
ming from that interaction on the background of which 
develops the fluctuating induced new attraction between 
the surfaces which is in fact the critical Casimir force. 



In the definition (1.1) used here, that is the common one 
when one considers short-range systems, these both effects 
are superposed simultaneously. Therefore, here, generally 
speaking, one should expect a crossover from a regime 
governed by the critical Casimir force (in the sense of a 
f luct uation induced force; it is of the order of L~ d , see Eq. 
(L3)) to the one govern by the direct attraction (of the 
order of i _<T ; note that if d — a they will of the same 
order being dominating in different temperature regions). 
An interesting case when forces of similar origin are act- 
ing simultaneously is that one of the wetting when the 
wetting layer is nearly critical and intrudes between two 
noncritical phases if one takes into account the effect of 
long-range correlations and that one of the long-range van 
der Waals forces p2|,p3|]. 

In quantum systems additional new features will be 
observed since the "temporal direction" corresponds for- 
mally to a short-range type interaction in the correspond- 
ing classical analog of the system, i.e. one unavoidable has 
"anisotropy" in the spectrum of a quantum system with 
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long-range interactions. The effects depend on which di- 
mension - the temporal or the spatial one - is the finite 
one. When the spatial dimension is finite we only mention 
here and will demonstrate in the current article, that an 
effect similar to that for the classical systems exists. If the 
finite dimension is the temporal one such effect will not 
be observed since then the long-range interaction is inter- 
surficial and is not effected directly by the "finite size" of 
the system. 

The plan of the paper is as follows. In Section 2 we 
define some generalization of Casimir amplitudes in quan- 
tum systems and present some hypotheses for the corre- 
sponding excess free energies. Then in Section 3 we give 
a brief review of the model and the basic equations for 
the free energy and the quantum spherical field in the 
case of periodic boundary conditions. Since we make use 
of the ideas of the FSS theory, the bulk system in the low- 
temperature region is considered as an effective (d + z) 
dimensional classical system with z finite (temporal) di- 
mensions. This is done to make possible a comparison 
with other results based on the spherical type approxi- 
mation, e.g., in the framework of the spherical model and 
the QNLcM in the limit n — > oo. The scaling forms for 
the excess free energy, the spherical field equation and the 
Casimir force are derived for a i<r < d < 1<t dimensional 
system with a film geometry in Section [I| In Section ^| 
we present some results for the Casimir amplitudes in the 
case of short-range interactions and in some special cases 
of long-range interactions. The paper closes with conclud- 
ing remarks given in Section 0. 



2 Casimir amplitudes in critical quantum 
systems 

Let us consider a quantum system with a film geometry 
L x co d_1 x L T , where L T ~ ^(fc^T) is the "finite-size" 
in the temporal (imaginary time) direction and let us sup- 
pose that periodic boundary conditions are imposed across 
the finite space dimensionality L (in the remainder we will 
set h = ks = 1). Let f(T,g,H;L\d) be the free energy 
density of this system. Then, according to the dimensional 
crossover rule, the Privman - Fi sher hypothesis pif for fi- 
nite classical systems and Eq. (1.2) in the quantum case 
one could state that 



lr(T,g,H;L\d) = (TL T ) L^ d +^ X^(x l7 x 2 , p\d), 



L 



with scaling variables 



5i = L^^Sg, x 2 = hL A ' u and p = L Z /L T . 



(2.1) 



(2.2) 



Here A and v are the usual critical exponents of the bulk 
model, h is a properly normalized external magnetic field 
H , Sg ~ g — g c , and X" x is the universal scaling func tion 
of the excess free energy. According to the definition (1.1) 
of the Casimir force, one obtains immediately 



Casimir 



(T,g,H;L) = (TL r ) L-^X^ BBimil (x u x 2 , p\d), 

(2.3) 



where the universal scaling functions of the Casimir force 
XQ asimir (xi, x 2 , p\d) is related to the one of the excess free 
energy X£ x = X£ x (xi, x 2 ,p\d) by 



Casimir 



1 dX u 

(x!,x 2 ,p\d) = -{d + z)X^ x - -xi-^r 



A 8X^ 

"! x %~5 

v ax 2 



zp- 



dp 



(2.4) 



It follows from Eq. ( |2.3| ) that depending on the scaling 
variable p one can consider the general case of Casimir 
amplitudes 



^Casimir (H<0 = *Ca 



(0,0, p\d). 



(2.5) 



The classical amplitudes A a ^(d) (for ( a, b) = periodic 
boundary conditions) introduced by Eq. (|l.3j ) are particu- 
lar cases of ^Casimir (p\d) f° r P = 0) i- e - T = (we remind 
that our system, according to the dimensional crossover 
rule, is formaly equivalent toarf+z dimensional classical 
one). 

In addition to the above "usual" excess free energy and 
Casimir amplitudes, denoted by the superscript "it", one 
can define, in a full analogy with what it has been done 
above, "temporal excess free energy density" / t ex , 

/ t ox (T, g, H\d) = f(T, g, H;oo\d) - /(0, g, H;oo\d) (2.6) 



and "temporal Casimir amplitudes" 
fr(T,g c ,0\d) = TL; d / z A 



Casimir 



(d) 



(2.7) 



Whereas the "usual" amplitudes characterize the lead- 
ing L corrections to the bulk free energy density at the 
critical point, the "temporal amplitudes" determine the 
leading temperature-dependent corrections to the ground 
state energy of an infinite system at its quantum critical 
point g c . 



If in Eq. (2.7) the quantum parameter g is in the vicin- 
ity of g Cl then one expects 

fr(T,g,H) = TL- d / z Xl (x\,x\\d) , (2.8) 

i.e. instead of the amplitude ^casimir (*0 one nas a scaling 
function X^ (x\,x 2 \d) which is the corresponding analog 
of X£ x (xi, x 2 , p\d). The scaling variables now are 



l}l vz 8g and x\ = hL A ' vz . 



Obviously 



Casimir 



(d)=*2*(o,o|d). 



(2.9) 



(2.10) 



Let us finally note that if z = 1 the temporal excess free 
energy introduced here coincides, up to a (negative) nor- 
malization factor, with the proposed by Neto and Fradkin 
p5| definition of the non-zero temperature generalization 
of the C-function of Zamolodchikov. 

Now we pass to study the quantities introduced above 
on the example of one exactly solvable model. 
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3 The model 

The model we consider is described by the Hamiltonian Pj 



H 



\d E v i - \ E J u>SiS v + l^J2 s e- H J2 S * 



where Se are spin operators at site I. The operators Vg 
play the role of "conjugated" momenta (i.e. [<S>£,6>^'] = 
0, [P t ,Pf] = 0, and [Pt,S t >] = iS U ', with fi = 1). The 
coupling constant g measures the strength of the quantum 
fluctuations (below it will be called quantum parameter), 
H is an ordering magnetic field, and the spherical field /i 
is introduced so as to ensure the constraint 



E<^ 2 > 



N. 



(3.2) 



Here N is the total number of quantum spins located at 
sites of a finite hypercubical lattice A of size L\ x x 
■ • • x Ld = N and (• ■ ■) denotes the standard thermody- 
namic average taken with the Ham iltonian H . In Ref. |fjj , 
the equivalence of the model (3T) and the quantum 0(n) 
nonlinear sigma model in its large n-limit is shown. 

Let us note that in the last few years an increasing 
interest in the spherical approximation (or large n-limit), 
generating tractable models of quantum critical phenom- 
ena, has been observed 0,|§|2|,||,||,|o|. There are 
different possible ways of quantization of the spherical 
constraint. In general they lead to di ffer ent universality 
classes at the quantum critical point |l^,^,^,^8). The 
commutation relations for the operators Si and Vp to- 
gether with the kinetic term in the Hamiltonian ( |3.l|) do 
not describe quantum Heisenberg-Dirac spins but quan- 
tum rotors as is pointed out in Ref. JlT| . Since the quan- 
tum rotors model has been widely exploited in the field of 
high-temperature superconductivity (see, e.g. |jl4| and ref- 
erences therein) we hope that the treatment of the model 
(3.1) presented below might be of some interest to those 
problems. 

For nearest neighbour interaction different low tem- 
perature regimes and finite-size scaling properties of the 
model are investigated in Ref. |^9| . 

The free energy of the model in a finite region A un- 
der periodic boundary conditions applied across the finite 
dimensions has the form BQl 



0f A (0,g,H) 




70- 



2w 2 (0; M )' 



2 sinh 



-W 



i/3w (g; /i) 



(3.3) 



Here the vector q has the components j 2 ^" 1 , • ■ ■ , 2 ™ d j , 

| — Lj 2 1 , ■ • ■ , L ' 2 1 | for Lj odd integers, and (3 is the 

1). In (3.3) the spectrum 
JM CT , < a < 2. 



2 ' '2 

inverse temperature (with ks 
is lj 2 (q; (i) = g(n + U{q)) with U(q) 
In the above expressions U (q) is the Fourier transform of 
the interaction matrix where the energy scale has been 



fixed so that £7(0) = 0. The supremum in Eq. ( |3.3D is 
attai ned at the solutions of the mean-spherical constraint, 
Eq. (U), that reads 



1 = 



t 



OO 



N E E 

rrt — — oo q 



<p + U(q)/J + b 2 m 2 (j) 



2 • 



(3.4) 



where we have introduced the notations: b — (2irt) /A, A = 
\J gl J is the normalized quantum parameter, t = T/J - 
the normalized temperature, h = H/J - the normalized 
magn etic field, and = /i/J is the scaled spherical field. 
Eqs. ( |3.3| ) and (3.4) provide the basis for studying the 
critical behaviour of the model under consideration. 

In the thermodynamic limit it has been shown ]l7t that 
for d > a the long-range order exists at finite temperatures 
up to a given critical temperature t c (X). Here we shall 
consider the low-temperature region for < d < |<r. 
We remind that and |tr are the lower and the upper 
critical dimensions, respectively, for the quantum critical 
point of the considered system. 



4 Scaling form of the excess free energy and 
the Casimir force at low temperatures 

For a system with a film geometry L x oo d_1 x L T (where 
icr < d < |ff) , after taking the limits L2 — > 00, • • • , Ld — » 
00 in Eq. ( |3.3j ) with L\ — L, we receive the following ex- 
pression for the full free energy density (see Appendix [A|) 



f(t,X,h;L\d,a)/J 



2(f> 



Xk r l j , „ 1 / 1 d x n 

\aL-( d +^ ~ [°° , 1 
- — 2^ / dxx a 2 



(Any 



_i Jo 



n=l 



n 

° XP 1 -~x 



xG» i_< 



v°l 2 L% 



A_fcd_^ / d 
a \/7r \ a 



°° Kd 



E 



00 oc 



(27T)" 



dz f z \ 



mrz 2 



x exp 



-zL° 



AzL c 



.(4.1) 



Here k d 1 = |(47r) 2 r(d/2), xd is the radius of the spher- 
icalizcd Brillouin zone, 



1 x 

G aJ3 (*) = 77= E 



r(k + 1/2) t k 

¥ r (ak + (3) fc! 



(4.2) 



(it was introduced in Ref. 

F„,a(y)= / x v+l J v {x)ex 9 {-yx' 7 )dx (4.3) 
Jo 
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and K v (x), and J u {x) are the MacDonald and Bessel func- 
tions, respectively. The main advantage of the above ex- 
pression, desp ite of its complicated form in comparison 
with Eq. (^3), is the simplified dependence on the size 
L which now enters only via the arguments of some func- 
tions. This gives us the possibility, as it is explained below, 
to obtain the scaling functions of the excess free energy 
and the Cas imi r force. 

In Eq. (|4.l| ) <j) is the solution of the corresponding 
spherical field equation that follow s by requiring the par- 
tial derivative of the r.h.s. of Eq. ( [4. 1| ) with respect to <\> 
to be zero. The bulk free energy /buik(i, A, h\d, a) results 
from f(t, A, h; L\d, a) by merely taking the limit L — > oo 
in it. Let us denote the solution of the corresponding bulk 
spherical field equation by (poo. Then for the excess free 



energy it is possible to obtain from Eqs. (1.2) and ( f4.1| ), in 
full accordance with Eq. (|2.l|), the finite size scaling form 



\r{t,\h-L\d,a) = (TL T )L-( d+2 )X c u x K,^,pM,a), 



L 



with scaling variables 

u _ Jl/V 



(4.4) 



L 1,u (X- 1 - X- 1 ) , x u 2 = hL A/v and p = L Z /L T , 



with L r 



X/t. Here the critical value of A = A c is 
1 



a; 



;( 27 



d d q(U(q)/JY 



(4.5) 

(4.6) 
-.a are 



and v^ 1 = d — ic, Ajv = i (d + , and z 
the critical exponents of the model jlTj. In Eq. (4.4) the 
universal scaling function X" x (x", x^Tpld, er) of the excess 
free energy has the form 



X: x (x u 1: x^p\d,a 
kd 



40TO- 

kd 



r 



yo) + Ux u 2 f(— 

2 \2/oc 



x l (Voo 

d 1\ / i+I 



1 

yo 



ySc 



r 



E 

rn—l 



(2y )" 



Kd, 



p J 



(22/oo) ' 



(J 



1 

4(4tt 



E 

n=l 



exp 



Ax J 



X%i-j {-x^yo) dx 



V2 



(27T) 



i±i E E 

n— lm— 1 



dz 



x exp 



-1,<T 

2 



(— ) 

Vm <T / 



4z,o 2 



By direct evaluation of the above expression it is easy 
to see that below A c if the finite system is in ordered phase 
(then j/o = J/00 = 0) X" x ~ that reflects the domi- 

nating contribution of the direct inter-spin long-rang e in - 
teraction in that region [Q. As we see from Eq. (4.4) 
one observes a crossover from L~( d+Z ^ behavior, where 
the fluctuation induced interactions dominate, to L~ a one 
where the direct inter-spin interactions become essential. 

For the Casimir force one obtains 



Fc'^(T,\h;L) = (TL T ) L^^^S™^ ,x% , p\d,a), 

(4.8) 

where the universal scaling function of the Casimir force 
-^Casimir ( x i i x 2 > pMj ct ) is related to that one of the excess 
free energy A c u x = X^, x$, p\d, a) by Eq. @. 

The above expressions for the scaling functions of the 
excess free energy (and the Casimir force) are the most 
general ones, which gives the possibility of a general anal- 
ysis including issues as: i) the sign of the Casimir force; 
ii) monotonicity of the Casimir force as a function of the 
temperature; Hi) the relation of the excess free energy scal- 
ing function to the corresponding finite-temperature C- 
function and its monotonicity properties; iv) finite-system 
generalization of the finite-temperature C-function, etc. 

In the present article we will concentrate on evaluation 
of the Casimir amplitudes for some special cases where one 
can obtain simple analytical expressions for them. 

It is clear that just due to the dimensional crossover 
rule L plays the same role for the finite system at t = 
as L T for the corresponding infinite quantum system. 
Therefore, by a symmetry that obviously arises when a = 
2, one should expect that the behavior of the two types of 
amplitudes ( "normal" and "temporal" ) will be essentially 
the same. We will see that explicitly below. 

For the model we study here one can show that 



Xl,{x\,x\\d,a) = \x\ (2/00 - 2/0) + \{x\f ( — - — 
2 2 \2/oo 2/0 



d 1 



(d/z + l)/2 _ (d/ Z + l)/2 

hfOC 



*V(£W*i£; liiit^L. (4 . 9) 



d r 1 



(4 



(see Appendix ^ for details of the calculations), where 
j/o = 4>L a and j/oo = 4>ooL a . 



m=l [myfyo 
Here the scaling variables are defined by 

x\ = L X J V ' (~ - i J , x\ = hL?/*», (4.10) 

2/0 = L^tfio and j/oo = -^r^oo- Note that j/o is the solution of 
the corresponding spherical filed equation for the nonzero- 
temperature system, whereas 2/00 is the solution for the 
zero-temperature (" infinite" in t he "temporal" dimension) 
one. The direct evaluation of (19), supposing the finite 
system in a ordered state, shows that X^ x remains of the 
^)same order in that region. This is exactly the behavior to 
be expected, as we mentioned in the Introduction, despite 
the long-range nature of the interactions, since the finite 
"dimension" of the system is now the temporal one. 
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Now we are ready to investigate in a bit more detail 
the behavior of the Casimir amplitudes as a function of 
d, a and p. 



5 Evaluation of Casimir amplitudes 

In this section we determine the Casimir amplitudes of 
the model in the case of short-range interactions at d — 2 
and in the special case d = a of long-range interactions. 



5.1.2 One-dimensional system with long-range interactions 

(d = a=l) 

We can o btai n a relatively simple analytical expression 
from Eq. (4.7) only in the particular case a = 1. In this 



case the functions G and T become 

Gi/ 2 ,3/4(-z) = ^ exp (z 2 /2) K 1/4 [ Z -\ ( r,.r, ! 



and 



5.1 "Usual" Casimir amplitudes 

5.1.1 Two-dimensional system with short-range interactions 

(d = a = 2) 

In this case essential simplifications in the expression for 
the Casimir forces can be made. The functions G a ,p and 
^a,0 used in the general expression of the free energy (4.1) 
turn into the explicit forms 



1 



and 



Gi,i/ 2 (z) = ^exp 
V 71 " 



F a A z ) = ( 2z ) a 1 ex P 



4z 



(5.1) 



(5.2) 



At the quantum critical point A = A c , h — this leads to 

( yoo = o ) 



JC»(Q,0,p|2,2) 



1 3/2 



12tt 



(5.3) 



where the primed summation over the integers m and n 
indicates that the term corresponding to m — n = is 
excluded. 

Since it is not clear how to obtain an explicit analyti- 
cal solution for yo i n a fil m geometry at nonzero temper- 
ature, the above expression cannot be simplified further, 
but has to be analyzed numerically (see, e.g., |29| for a 
numerical analysis of the spherical filed equation^Never- 
theless, the above expression can be significantly simpli- 
fied at zero temperature. Then one shows that the solu- 
tion yo of the spherical field equation for the finite system 
with a film geometry L x oo x L T , at zero temperature 
(i.e. 1 -C L <C oo, L T = oo) and at the quantum critical 
point A = A c , h = gMis y = 4 In 2 (V5/2 + 1/2). 
Setting this value of yo in ( |5.3| ), taking into account that 
K3/2( x ) — \J t^I (2ic) exp(— x)(l + l/x), and using the prop- 
erties of the polylogarithm functions Lip(x) | p3[ ], we obtain 
after some algebra that the Casimir amplitude is 



Z\£ asimir (0|2,2) 



2C(3) 
5tt 



-0.1530. 



(5.4) 



Here C{ x ) i s the Riemann zeta function. 



T u ,x{y) - ^r(^ + 3/2) 



(l + y2) 



I/+3/2 : 



(5.6) 



respectively. In order to obtain explicit results for the am- 
plitudes, numerical evaluations are unavoidable even in 
the simplest case correspondi ng to zero temperature. In 
this case we obtain from Eq. (|4.7|) 



AC(0,0,0|1,1) 



. OO 

s/yo 

°7T 2 <_i JO 



dxx 2 exp 



Ax 



xKi 



3tt 



y 



3/2 



(5.7) 



Here yo is the solution of the equation for the spheri- 
cal field which can be obt aine d by requiring the partial 
derivative of the r.h.s of ( |5.7|) with respect to yo to be 
zero. Solving the last equation numerically we end up with 
yo = 0.6248. After substitution of the solution in Eq. (5.7) 
we obtain for the Casimir amplitude 



A 



Casimir 



(0|1,1) = -0.3157. 



(5.8) 



The result given by Eq. (5.8) shows that the Casimir am- 
plitude in the case a = 1 is of larger magnitude than the 
one in the case of short-range interaction. One can ask 
whether this is just a coincidence or the Casimir ampli- 
tudes are increasing function of a for a given fixed d. 



5.1.3 Relation with the Zamolochikov's C-function 



In terms of the critical-point value of an analog of the 
finit e-temperature C-function p5| , the result given by Eq. 
(5.4) can be rewritten in the form 



^ asimir (0|d,2) 



-n u (d, 2)c u (d,2), 



(5.9) 



where in analogy with the short-range interaction case 
above we define a number c w via the relation 



c u (d,a) 



.An / 
Casimir ' 



\d, a)/n u (d, a) 



(5.10) 



Here, as usual, the normalization factor n u (d, a) is chosen 
so that to ensure c a = 1 for the corresponding Gaussian 
model, i.e. 



n u {d,a) 



2 g/2 ^(| + f)C(rf- 
7r<Z/2|r(-f)| 



(5.11) 
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For (j = 2 we immediately obtain c"(2,2) = 4/5, where 
n u (d,2) = r({d + l)/2)C(d + 1)/tt^ +1 ^ 2 becomes the 
normalization factor given in [p5[| . In that way we repro- 
duce the well known result for c due to Sachdev Q who 
considered an example of a three dimensional conformal 
field theory. This coincidence of the values of c is due to 
the fact that both models belong to the same universal- 
ity class. For more details on the behavior of the finite 
temperature C-function in the case a = 2,d = 1,2,4 
see, e.g. Q. The d-dependence of the value c u (d, 2) = 
z4£ asimir (Ojd, 2) /n u (d,2) has been considered in |?5) for 
d-dimensional (2 < d < 4) conformally invariant field the- 
ory. The relation between the C function and the Casimir 
force for the classical version of the model and short-range 
interaction has been analyzed in some details for d = 2 
in@. 

In the particular case of lon g-ran ge interact ion d = 
a = 1, in accordance with Eq. ( 5.10 ) and Eq. ( |5.8| ) one 

gets 

c u (l,l) = 0.606. (5.12) 



5.2 "Temporal" Casimir amplitudes 

5.2.1 cr-dimensional system (d = a) 

Let us note that unlike the case of " usual" Casimir ampli- 
tudes here it is possible to consider the more general case 
d = a, where < a < 2. From Eq. ( |4.S| ) one obtains a gen- 
eral expression for the temporal Casimir amplitudes for a 
system with geometry oo d x L T , at the quantum critical 
point A = A c , h — 0, 



^Casimir a ) 



+i 



kd r f d 
aJU \ a 



E 



Kc 



(myo) 



(myo) 



d [l 



(5.13) 



where the scaling variable yo is the solution of the corre- 
sponding equation for the spherical field. We no tice here 
that in the particular case d/a = 1 Eq. ( pUg] ) simpli- 
fies considerably. That is why we are going to investigate 
namely this case. In a way similar to that explained in the 
case of short-range interactions, one obtains (0 < a < 2) 



^Casimir ( a > a ) 



16 C(3) 



1 



5ct (4tt) ct / 2 r{a/2)' 



(5.14) 



Note that the defined "temporal Casimir amplitude" 
(a, a) reduces for a = 2 to the " no rma l" Casimir 
amplitude Z\£ asimir (0|2, 2), given by Eq. fTJ). This re- 



Casimir 



fleets the existence of a special symmetry in that case 
between the "temporal" and the space dimensionalities of 
the system. 

When a ^ 2 it is easy to verify that the following 
general relation 



^Casimir (°S °) 

M 

Casimir 



Stt 



(2,2) <j{A-Ky/ 2 r{<j/2) 



(5.15) 



between the temporal amplitudes holds. The r.h.s. of ( 5.15 ) 
is a decreasing function of a. 



5.2.2 Relation with the Zamolochikov's C-function 

As it has been already mentioned in S ectio n 2 if z = 1 the 
temporal excess free energy (see Eq. (2^)) coincides, up 
to a negative normalization constant, with the nonzero- 
temperature generalization of the C-function of Zamolod- 
chikov proposed in p5[ . For z / 1 a straightforward gen- 
eralization of this definition can be proposed at least in 
the case of long-range power-low decaying interaction 



C(T,g\d,z) = _T-< 1+d /^ 



v d/z 

n{d, z) 



fUT,g), (5.16) 



where z = a/2, the nonuniversal constant v in our nota- 
tions is v = TL T (see Eq. ( |2.8[ )), and the normalization 
factor is taken to be such that the corresponding Gaus- 
sian model with the considered type of interaction will 
have C(T, g c \d, z) = c t (d,a) = 1 at its critical point, i.e. 



n*(d,a) = - C(1 
a 



2d/a) r(2d/a) 
(47r) d / 2 r(d/2) ' 



(5.17) 



Let us note that the above choice of the normalization 
constant n*(d,z) preserves not only the c value for the 
Gaussian model, but also the corresponding one fo r the 
spherical model if d = a. Indeed, in that case from ( |5.14 ) 
and the above definitions one again obtains that 



c\a,a) = 4/5 



(5.18) 



for the spherical model. The result given by Eq. ( 5.18] ) is a 
generalization to the case of long-range interaction of the 
Sachdev's result. 



6 Concluding remarks 

In the present article the free energy of a system with 



a geometry L x ex 
derived (see Eq 



d-l 



L T (where \a < d < |cr) 



is 

( |4.l| )). For a = 2 this general result re- 
duces to the one reported in J3(| where only the case of 
short -ran ge interactions has been considered. The expres- 
sion (4.4) represents actually the verification of the ana- 
log of the Privman-Fisher hypothesis |24[] for the finite- 
size scaling form of the free energy (formulated initially 
for classical systems) in the case when the quantum fluc- 
tuations are essential. Note, that in that case one has a 
finite space dimension and one additional finite dimen- 
sion that is proportional to the inverse temperature, which 
provide different types of critical regimes and Casimir 
amplitudes. According to the finite-size scaling hypoth- 
esis [|4|[3(| one has to expect that the temperature multi- 
plying the universal scaling function will be with exponent 
p = 1 + d/z pq| , where the dynamic-critical exponent z 
reflects the anisotropic scaling between space and "tem- 
poral" ( " imagina ry-time" ) directions. 



Eqs. (Hj7f|4 



present a general expression for the Casi- 
mir force in the quantum spherical model. In the classical 
limit (A = 0) for a system with short-range interaction it 
coincides with the corresponding one derived in ^,0] for 
the classical spherical model. 
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In order to derive the Casimir amplitudes in a simple 
analytical closed form, some particular cases {d = a) have 
been considered: 

1) For the short-range c ase {a = 2) the corresponding 
amplitude is given in Eq. (5.4). This amplitude is equal 
to the "temporal Casimir amplitude" for the 0(n) sigma 
model in the limit n — > oo ]33[ . In the short-range case we 
have demonstrated explicitly that the two models, due to 
the fact that they belong to the same universality class, 
indeed possess equal Casimir amplitudes as it is to be 
expected on the basis of the hypothesis of universality. 

2) In the long-range case (d = o ^ 2), the correction 
to the ground-state energy of the bulk sys tem d ue to the 
nonzero temperature is determined by Eq. (5.14). One ob- 



serves that in this case c (a, a) — 4/5 does not depend on 
a. This could be understood by noting that by changing 
a one does not change the exponent in the spectrum that 
correspo nds t o the "temporal" (finite-size) dimensionality 
(see Eq. Q). 

3) At zero temperature we evaluated numerically the 
Casimir amplitude for the particular cas e d = a = 1 of 
the long-range interaction. The result (|]^) shows that the 
Casimir amplitude is of a larger magnitude than in the 
case of short range interaction [a = 2). Furthermore, the 
universal amplitude c"(l, 1) is no longer c-independent, 
because the finite-size part of the spectrum is c-dependent 
in this case. 

In accordance with the general expectations, all the 
amp lit ude s t hat we hav e de rived are negative (see Eqs. 
(fj), fTJ), |[l|)) and |[|). 

Fi nall y, let us note that the basic expression (see 
Eqs. (4.7)) for the scaling function of the excess free energy 
can be used as a starting point for generalization of some 
of the existing results on the C-function to the c ase o f 
long-range interactions. We have suggested in Eq. ( 5.16 ) 
a generalization of the nonzero-temperature C-function, 
proposed by Neto and Fradkin in ]25| , to the case of power- 
law long-range interactions. For the quantum spherical 
model this definition leads to c = 4/5 for any d = a which 
generalizes to long-range interactions the corresponding 
result for the case of short-range interactions (d = a = 2) 
due to Sachdev [B3|. 
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A Mathematical Appendix 



First we explain how from Eq. ( p.3[) one can obtain (4.1) 
for a syste m wi th a geometry L x oo d_1 x L T . It is easy 
to see that (3.3) can be rewritten in the form 



f(t,X,h;L)/J 



t 

N 



^2 m 1 2 smh 

q I 



h 2 4> 

24>~ 2 



A(L)-tJ2u n (L), (A.l) 



where 



and 



exp 



-n-V^+|qF 



(A.2) 



(A.3) 



In order to calculate A(L) we use the identity [ |3l[ 

VTT^ =- [l - exp (-zt)} r x - a i 2 G a ^ % {-n 

1 Jo 



-1 



(A.4) 



(see Eq. (1^2) for the definition of G a .p). Since we are 
interested in the geometry L x oo d ^ 1 x L T , one has by 
standard arguments 

i i r l lL ~ 1]/2 

JV q ( 27F ) J %=-[L-l]/2 

In the last sum over q, by using the asymptotic formula 



exp 



- a {-) 




>(A.6) 



1=1 



valid for L >> 1, we obtain, after some algebra, 

A(L) = A(oo) + SA(L), (A.7) 



where 
A(oo) 

and |l| 

SA(L) = 



dqi 



" dq d \f<t>+(q! + ... + q 2 d r /2 

-7T 

(A.S 



A a x ^ 

7,, xd/2 2-^1 



(4tt 



1=1 



x -a /4-d/2-l exp 



l 2 L 2 
Ax 

(A.9) 



Since the only singularities of A(oo) as a function of 4> are 
coming from small g's, it is justified to use a sphericaliza- 
tion of the Brillouin zone, which leads to 



ai ^ x , f XD dx n. — 



Xk 



d „d 



2d 



x D {xl 



jl/2 

d x°, 



(A.10) 
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where 2-F1 is the hypergeometric function. Now it is clear 
how the "first half" of ( }4.l| ) can be obtained. Next we 
t urn to evaluation of the term U n {L). Taking into account 
( A. 5 ), we rewrite U n {L) in the form 



U n (L) 



L -i (i-i)/2 k 
TTT^T E / 



q=-(L-l)/2 



x exp 



-nj\Jcf> + (qt + ... + q 2 d y 



dqd 
(A.12) 



Using the Poisson summation formula 

b 00 „b 

/ (n) = / dn exp [i27rfcn] 

n=a k=-oo a 

+\[f(a) + f(b)\ 
we obtain from the above expression 

U n (L) = U n (oo) + 8U n (L), 

where 

A 



f(n) 



U n (oo) 



dxx d 1 exp 



n—y <j) + x° 



(A.13) 



(A.14) 



(A.15) 



and 



SU n (L) 



~7d E / d<l1 ■ ■ ■ / dqd C ° S ^ llL ^ 

7r) j =1 J— 7T J — 7T 



(2 

x exp 



-n-y l (f ) +(q 2 + ...+ 



A.16) 



I n the low-temperature limit t « 1 one can replace in 
(A.15) xd by infinity. Then, using the integral represen- 
tation of K u (x) 

K u (2y/zi) = K- v {2-fzi) 

fOO 



1 



we derive 
U n (oo) 



A kg ( d 



te - - ) dx(A.17) 

xJ 



d + l I A , 1 



2t ' 



(A.18) 



We are left to deal no w onl y with SU n (L). Sphericalizing 
the Brillouin zone in ( A.16| ), performing the integrations, 
by using the integral representation for the Bessel function 

Ju{z) 



(a 2 — x 2 ) 13 1 exp [iXx] dx 



v~n ( y 



18-I/2 



J /3 _i /2 (aA)(A.19) 



( Re > 0) we get 

2L -d/2+l 



7TTd72- £ / TdT^T J d/2-i (^) 
(2tt) i=1 Jo ' 



x exp 



—n—y 1 <p + x° 



(A.20) 



In the low-temperature limit the upper limit of integration 
i n the above expre ssions can be replaced by infinity. From 
( A.18| ) and (A.20) one obtains the last two terms in Eq. 
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